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The arrest of Langmuir-wave collapse by quantum effects, first addressed by Haas and Shukla 
[Phys. Rev. E 79, 066402 (2009)] using a Rayleigh-Ritz trial-function method is revisited, using 
rigorous estimates and systematic asymptotic expansions. The absence of blow up for the so-called 
quantum Zakharov equations is proved in two and three dimensions, whatever the strength of the 
quantum effects. The time-periodic behavior of the solution for initial conditions slightly in excess 
of the singularity threshold for the classical problem is established for various settings in two space 
dimensions. The difficulty of developing a consistent perturbative approach in three dimensions is 
0^ ' also discussed, and a semi-phenomenological model is suggested for this case. 

o i 

O ! PACS numbers: 52.35.Mw, 52.35.g, 52.65.Vv 

(N 



> 
O 



- 1—1 

X 



I. INTRODUCTION 



o 

Special interest was recently devoted to quantum corrections to the Zakharov equations for Langmuir waves in a 
plasma First considered in one space dimension Qjthe model was then extended to two and three dimensions 
p ^ ' in a formulation retaining magnetic field fluctuations [4J . In a non-dimensional form, the equations that govern the 
' amplitude E of the electric field oscillations and the number density n read 

r^* I idtE - aV x (V x E) + V(V ■ E) 

= nE + TVA(V-E) (1) 
d tt n-An = A|E| 2 -TAV (2) 



In the above equations, the parameter a defined as the square ratio of the light speed and the electron Fermi velocity 
is usually large. The corresponding term is nevertheless moderate, because, magnetic effects are relatively weak, 
making E close to a gradient field. In contrast, the coefficient T that measures the influence of quantum effects is 
usually very small. We refer to [|| for a discussion of the physical regimes described by the present model and an 



00 ■ estimate of the plasma parameters. Typically, for a hydrogen plasma, one has a w 8.HF n 2/3 and T » 6.I0 6 n 1/3 , 
leading to a w 8.10 2 and T ss I0 -4 in the case of the rather high equilibrium densities no = 10 32 m~ 3 . 
Equation ^ originates from the hydrodynamic system 

8 t n + V • v = (3) 
O' <9 t v = -V(n + |E| 2 ) + TVAn, (4) 



governing the ion sound waves. 

For r = 0, it was shown [5|, [6|, [7|, [8| that for a ^ 1 and "small enough" initial conditions, the solution remains smooth 
for all time. In two dimensions, the smallness condition reads |Eo|l22 < |-R| 2 2 ~ 1-86 (where R is the ground state 
defined in (I25p ) and is optimal. In three dimensions, it requires that the plasmon number M and the Hamiltonian 
H (defined in ©,([Tn])) satisfy Af\H\ < 2.6 1(T 4 together with |VE || < \H\, conditions that are probably much too 
strict. 

Although not rigorously proved, the phenomenon of wave collapse is expected for T = 0, when the initial conditions 
are large enough. The existence of a finite-time blow up is indeed suspected when the Hamiltonian is negative, on 
the basis of numerical simulations and heuristic arguments (see Q for review). 

The question then arises of the possible arrest of collapse by quantum corrections. This issue was addressed in [3j by 
implementing an approach based on the Rayleigh-Ritz trial function method, in regimes where the quantum Zakharov 
equations can be reduced to a vector nonlinear Schrodinger equation. The latter results from the assumption that 
the density is slaved to the electric field oscillations (adiabatic approximation) 

- An = A|E| 2 -TA 2 n (5) 

and, because of the smallness of T, can be expressed to leading order as n = — |E| 2 - TA|E| 2 . Although this approach 
led to interesting conclusions such as the arrest of collapse by quantum effects and its replacement by a time-periodic 
solution, it nevertheless involves possibly questionable assumptions. For small enough T, quantum effects become 
relevant only very close to the singularity when the adiabatic approximation, even if valid at early times, hardly 
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holds. The rates of blow-up of the solutions of the Zakharov equations (with T = 0) are indeed such that all the 
terms in cq. ^ have the same magnitude in two dimensions, while <9 tt n ^> An in three dimensions (supersonic 
collapse). Furthermore, the Rayleigh-Ritz trial function method used to reduce the problem to the evolution of a 
few scaling coefficients, is based on an arbitrary choice the functional form of the solution. The aim of the present 
approach is to revisit the issue of collapse arrest by quantum corrections, mainly in regimes amenable to a systematic 
analysis. The paper is organized as follows. Section II provides a rigorous proof of the arrest of Langmuir collapse by 
arbitrarily small quantum effects, in the general framework of the Zakharov equations (|1I2|) . both in two and three 
space dimensions. Section III reviews the electrostatic approximation that is valid when the plasma is not too hot, 
as well as the so called scalar model that in the case of rotational symmetry does not prescribe a zero electric 
field at the symmetry center and a ring (2D) or shell (3D) profile for the electric field intensity. Section IV provides 
an asymptotic analysis of the dynamics in the presence of weak quantum effects for the scalar model in two space 
dimensions. This analysis is extended to the two-dimensional electrostatic equations with rotational symmetry in 
Section V. The difficulties of the three-dimensional problem are discussed in Section VI, where a phenomenological 
model based on a heuristic extension of perturbative calculations is presented. Section VII briefly summarizes our 
conclusions. 



II. ARREST OF COLLAPSE BY QUANTUM EFFECTS 



In this section, we present a rigorous proof of the absence of wave collapse for the general Zakharov equations with 
quantum effects (TTJ [5]) in space dimension d — 2 and 3. For this purpose, we first define the L p (R d ) and Sobolev 
spaces H s (R d ) as the spaces of scalar, vector or tensor functions respectively equipped with the norms [ll| 

\f\ LP = (j \f\ p d X ) 1/P (6) 

ijv = ( /(i+fc 2 n/(fc)iW) 1/2 , (?) 



where f(k) denotes the spatial Fourier transform of the function /. We will also use of a special case of the Gagliardo- 
Nirenberg inequality (see e.g. [l2j for the general formula), in the form: for / G {L q (R d ), V/ G L 2 (M. d )}, one has 

l/k»<^|v/|£ 2 |/|i; 9 , (8) 

where if is a positive constant whose optimal value is given in [l3| in two space dimensions and in a more general 

setting in [14]. In dimensions d = 1 and 2, 1 < q < p, while for d > 2, 1 < q < p < — — -. In all the cases, 

_ 2d(p - q) 
p[2d-q(d-2)]' 

Among the conserved quantities of the quantum Zakharov equations (|H2p , the number of plasmons and the Hamil- 
tonian play an important role in the regularity properties of the solution. They read [3j 

Af= J |E| 2 dx (9) 

1 2 1 

2" + 2 

+n|E| 2 + r|V(V • E)| 2 + | |Vn| 2 }dx. (10) 

Among all the terms arising in the Hamiltonian, only J n|E| 2 <ix could be non-positive and thus needs to be 
estimated. One has (denoting by C different constants) 



H = I \ IV • El 2 + a|V x El 2 + -n 2 + -Ivl 2 



J n |E| 2 dx ^ |n| L 4 ||E| 2 | L4/3 = \n\ L , |E| 2 8/3 



< C\n\ H i |E| 2 s/3 < \\n\ 2 H1 + £|E|i. /3 , (11) 

where the successive inequalities result from the Holder, Gagliardo-Nirenberg and Young inequalities. Using again 
the Gagliardo-Nirenberg inequality, we write 

|E| LS/3 s;C|VE|^ 2 8 |E|^ d/8 . (12) 
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It follows that 



nlEl^x^M^+^-^IVE^ 2 . (13) 



It is then convenient to rewrite the Hamiltonian in the form 



H = | VE|| 2 + (a - 1)|V x E|| a + r|V(V • E)| 2 L 

r 

2 



r, -| 2 ffl + (^-^)H 2 L 2 + i|v|| 2 +/n|E| 2 dx. 



Using eq. (|T3|) . one gets the upper bound 

|VE|| 2 + (a - 1)|V x E| 2 L2 + T|V(V • E)|£ 2 



(14) 



+ \\n\ 2 m + (^-\)\n\l. + -\v ]L2 

^iWl + Jlnl^ + ^-^IVEl^ 2 . (15) 

Since a and T are respectively larger and smaller than 1, one in particular has 

|VE|| 2 < \H\ + ^M 2 - d/4 (\VE\ 2 L2 ) d / 4 . (16) 

which implies that, for any T > 0, | VE|^ 2 remains uniformly bounded in time. Equation (|15[) also provides an uniform 
bound for 7T,| jj-i , |v|^ 2 and |V(V • E)| 2 2 . One then easily derives, using standard methods, the existence for all time 
of a classical solution for the quantum Zakharov equations, both in two and three dimensions. 



III. ELECTROSTATIC LIMIT AND SCALAR MODEL 



The large value of the parameter a makes the magnetic fluctuations actually subdominant, leading to the so called 
electrostatic approximation [l|. For this purpose, it is convenient to derive from eq. (|TJ) the system 

id t {V ■ E) + A(V • E) = n(V • E) + V?i • E 

+rA 2 (V-E) (17) 
id t (V x E) + aA(V x E) = n(V x E) 

+Vn x E. (18) 



Even if the initial electric field is a gradient, V x E is driven by the last term in the right hand side of eq. (Q2 
Nevertheless, when a is large, a stationary-phase argument applied to this equation, easily shows that V x E saturates 
at a level that scales like 1/a. Thus, although small, it contributes to eq. ([T]) but not to the Hamiltonian that, as 
a — > oo, has a finite limit obtained by neglecting the term involving the coefficient a. 

Writing E = — V^> H — Ei and substituting in eq. (fTTl) . one gets to leading order 
a 

A(id t 4> + AVO = V ■ (nVt/j) + rA 3 V>, (19) 
dun — An = AQVV'I 2 ) - TA 2 n. (20) 

A rigorous proof of the convergence is given in [l5| in the cases where the solution is globally smooth. The proximity 
of a singularity is nevertheless not expected to affect the ordering between the solenoidal and gradient components of 
the electric field. 

The analysis of the solution near collapse is often performed, assuming that the fluctuations involve rotational 
symmetry [lj. In this case, introducing E = —dip/dr, eqs. (IT9|) - (l20l) reduce to 

idtE + A^E^nE + TAl^E (21) 
d tt n- A r n = A r \E\ 2 -TA 2 r n (22) 
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where Af^ = d r r^^ d ^ d r r d ~ l and A r = r~( d ~ 1 ^d r r d ~ 1 d r . These equations are supplemented with the boundary 
conditions E(0, t) = <9 r n(0, t) — E(oo, t) = n(oo, t) = 0. 

It was early noticed [10| that the assumption of an electric field vanishing at the center of symmetry (taken as the 
origin of coordinates) with a shell profile for the intensity profile is hardly consistent with a realistic model. Relaxing 
this assumption of zero electric field at the center while retaining an isotropic intensity profile of the electric field, is 
not possible when the detailed dynamics are retained. It may thus be suitable flG ] to abandon the vector character 
of the problem, in order to preserve a non-zero electric field at the center of the cavity as suggested by numerical 
simulations of the vector Zakharov equation near coll apse 11611 . while keeping the rotational symmetry necessary for 
implementing an asymptotic analysis in the spirit of [171 . 1 1 81 ] . We are thus led to consider the influence of quantum 
effects in the framework of the "scalar model" 10] 

id t E + AE = nE + TA 2 E, (23) 
d tt n -An= A\E\ 2 ~ TA 2 n. (24) 

with the condition that E and n vanish at infinity and satisfy d r E(0, t) = d r n(0, t) = 0. Isotropic solutions of these 
equations are not only stable but also attractive near collapse [19[ . It is interesting to notice that direct numerical 
simulations [l6| of the collapsing solutions of the vector Zakharov equations {T][2]) with T = indicate that the 
anisotropy is in general rather moderate and the rates of blow up identical to those of the scalar model fl9l ] . 

For r = and in the adiabatic limit where n = — |-E| 2 , isotropic solutions of eq. (|23p identifies with the vortex 
solutions corresponding to a rotational number m = 1 , of the (scalar) two-dimensional nonlinear Schrodinger equation 
[20| . Taking the adiabatic limit with r ^ 0, one gets a non local extension of the cubic nonlinear Schrodinger equation 
with both second and fourth order dispersions studied in [2l[ and [22| . 



IV. ASYMPTOTIC BEHAVIOR OF THE SCALAR MODEL IN TWO DIMENSIONS 

In the analysis presented in Q, the space dimension has no qualitative effect, but this is not the case in the 
framework of a systematic perturbative approach that, to be fully consistent, requires a small expansion parameter 
usually associated with closeness to a critical regime. 



A. The classical regime (T = 0) 



The two-dimensional regime deserves a special attention because it is amenable to a detailed mathematical analysis. 
It was indeed proved that for two-dimensional smooth initial conditions such that the initial density no £l 2 , dtfio S 
H^ 1 and the initial electric field E e H 1 obeys the condition |.Eo|.l 2 ^ \R\l 2 > where R is the unique positive solution 
(ground state) of 

AR - R + R 3 = 0, (25) 

the solution of the classical scalar model (T = 0) exists for all time in these spaces 0] and is unique [23|. Optimal 
local existence results in spaces of very weak regularity appear in [24j . Further regularity properties of the initial 
conditions are also preserved in time. The ground state R is radially symmetric [25| and obeys the relation J (| Vi?| 2 — 
(l/2)R 4 )rdr — 0, where the right hand side can be viewed as the Hamiltonian of the standing wave solution e lt i?(|x|) 
of the nonlinear Schrodinger (NLS) equation iip t + Atp + |V'| 2 ^' = 0. It is furthermore interesting to notice that, unlike 
the two-dimensional NLS equation, the Zakharov system does not have blowing up solutions of minimal L 2 -norm. 

Although there is no rigorous proof of existence of a finite-time singularity for larger initial conditions, one has 
the following result [7| (valid both in two and three dimensions). Suppose Ti < and that the solution (E,n,\r) is 
radially symmetric. Then either + \ii\l 2 + I v |l 2 —> 00 as t —> f+ with finite, or {E, n, v) exists for all time and 

+ \n\ L 2 + \~v\l2 — ► oo as t — > oo. Numerical simulations clearly indicate the stability of isotropic solutions and 
that, near collapse, every solution becomes locally isotropic pj| . Furthermore, in two dimensions, there exist exact 
self-similar solutions of the classical scalar model that blow up in a finite time, of the form 



E(r,t)= 1 p(- 
an i« -i v< 



»( e+- 

e 



g(t.-t) Mt,-t) 



(26) 

n(r '* ) = ^CT M b^)' (27) 
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under the condition that (P, M) satisfies the system of ordinary differential equations in the radial variable 

AP-P- MP = 0, (28) 
al(r} 2 M m + 6r)M v + 6M) - AM = AP 2 . (29) 

The parameter ao entering the self-similar solution is not universal. One has the following results [26[ for existence of 
solutions to (|2"5|) - (f2"9"]l . There exists > such that Va with < a < Oq , there is a solution (P ao ,M ao ) S H 1 x L 2 
of (gSD-dMD with P ao > 0. Furthermore, when a -> 0, (P ao ,M ao ) tends to (R, -R 2 ) in H 1 x L? and for all c > |i?| L a, 
there exists ao c > 0, such that for all ao with < ao < ao c there is a unique solution (P ao , M ao ) with P ao > and 

l^«oli»<C 

Numerical simulations of the classical scalar model were performed in [l9( where it is shown that for a smooth 
initial condition with |Po|l 2 > |P|l 2 j the solution of the initial value problem approaches the self-similar blowing-up 
solution. Furthermore, in a series of simulations with initial conditions such that |-Eo|l 2 approaches \R\l 2 from above, 
it was observed that the estimated value of the parameter ao monotonically decreases to zero (see table 1 ofpjjj]). This 
result is consistent with the existence of a sequence of blowing up solutions 



E » - T^' ( ^ +T4p) ^(T^i)> < 3 »» 
= (T^Mi^l) (31) 

obtained by choosing f* = —9 = a^ 1 in eqs. (|2"6")) - (f2"T)) . with the profiles obeying (j2"5)) - (f2"9")) and thus converging to 
(R, —R 2 ) as ao — > 0. For ao = 0, the corresponding solution is smooth, consistent with the regularity of the solutions 
such that \E\ L 2 = \R\ L 2 Q. 

The above observation suggests a pertubative analysis of the influence of quantum effects for initial conditions such 
that Po|l 2 is slightly above the threshold for collapse (27j . 

B. Influence of quantum effects 

The method for constructing a solution in the presence of weak quantum effects is to assume that these perturbations 
induce small corrections of the self-similar profile of the collapsing solution, but modify the scaling parameter A whose 
time evolution is prescribed by the conservation of the plasmon number and of the Hamiltonian. As noted in [l8| , an 
alternative variational approach based on the existence of a Lagrangian density is a priori possible, using the actual 
perturbative expansion of the fields. We shall not follow this direction here and concentrate on direct expansions. 

Using a dynamical rescaling transformation, we first define E(r, t) = A _1 J7(^, r), n(r, t) = A~ 2 iV(£, r), 
n t (r,t) = \~ 3 W(£,t) and v(r,t) = \~ 2 V((,,t), with £ = r/A and r = J* A(s) _2 ds. Introducing a = -A t = -A~ 2 A r , 
the rescaled quantities obey (A = ^~ 1 d^^d^ ) 

i[U T + a\(U + £Ut)} + AU — NU = T\- 2 A 2 U (32) 

N T + a\(2N + £Ns) = XW = -X^d^V) (33) 
V T + a\(2V + + XN 6 

= -Xd^Ul 2 + TX^d^AN (34) 
W T + aX(3W + £W{)- XAN 

= XA\U\ 2 -TX~ 1 A 2 N. (35) 

It is then convenient to write U = e lT e~ laX ^ and to define b — (aX) T + (aX) 2 . Equation (|32p is replaced by 

iS T - S + AS - NS + b^-S = — e iaA « 2 / 4 A 2 (e- wA « 2/4 5). (36) 

4 A 2 

The assumption is then that we are looking for solutions whose time dependency is only through the scaling parameter 
A, and thus also through the functions a and 6. This leads to neglect the contributions S T , N Tl V T and W T in the 
above equations. Combining the resulting equations for N and W, one gets 

a 2 £(AH - AN = AIS1 2 - ^A 2 iV. (37) 

X z 
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FIG. 1: Profiles of the positive solution (ground state) of the equation AR — R + R 3 = and of the density N = — R 2 . 



where the operator C is defined by C — £, 2 d^ + 6£c?£ + 6. 

For r = 0, the self-similar solution is reached as t approaches the singularity time i* (i.e. r — > oo) and a tends to a 
constant ao, which leads to the conclusion that in this regime A(t) = ao(t* — t). Furthermore, as mentioned in Section 
IV A, ao tends to zero when considering a sequence of initial conditions for which I-E0H2 approaches |-R|?2 ~ 1-86 
from above. In this limit, S is close to R, and N close to — R 2 . 

We now consider the effect of a small but non zero T, while keeping the assumption that the initial conditions are 
such that |23o|^2 is slightly in excess of |i2|| 2 . In this regime, the quantities a and b are supposed to remain small but 
are allowed to be time dependent. Furthermore, since the profiles of the various fields remain close to those of the 
blowing-up regime, we are led to expand them in terms of the small parameters a, b and T/X 2 in the form 

p 

S = R + a 2 <Ji + ba 2 + ^3 H (38) 

N = -R 2 + a 2 v l + bv 2 + ^-v i + --- (39) 
V = av x + --- (40) 

For sake of completeness, the leading order profiles R and — R 2 are plotted in Fig [TJ 

At the center of symmetry, d^ai, d^yi and vanish. All the functions also decay at infinity. This leads to the 
sequence of equations 

Acri - <7i + R 2 ai -Rui = (41) 

-Avi - 2A(Ra 1 ) = £(R 2 ) (42) 

C 1 d ( ($v 1 )=2R 2 + Z(R\ (43) 

£ 2 

Act 2 - <7 2 + R 2 a 2 - Rv 2 =~R (44) 
-Az/ 2 - 2A(i?cr 2 ) = (45) 

Acr 3 - a 3 + R 2 <7 3 - Rv 3 = A 2 R (46) 
-Aza, - 2A(i?CT 3 ) = A 2 (R 2 ) (47) 

Since the kernel of the operator A — 1 + 3R 2 is reduced to the null function under a radial symmetry assumption, 
all the above equations are solvable. We now write that the time evolution of the fields through the variations of the 
functions a, b and A are constrained by the conservation of the plasmon number and of the Hamiltonian that can 
both be estimated within the above perturbative expansion [l7|. 
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For the plasmon number, one has (up to a 2tt angular integration factor that we systematically omit) 

7V = y \E\ 2 rdr = 

J (it 2 + 2a 2 Rax + 2bRcr 2 + fdf + • ■ • . 



(48) 



The Hamiltonian reads 



H = J (\VE\ 2 +n\E\ 2 + l -n 2 + ^V 2 

+T\AE\ 2 + ^|Vn| 2 )rdr. (49) 

It involves |V-E| 2 = A [|ViS| + (a 2 A 2 /4)£ 2 S 2 ]. Substituting the expansion of the various fields and using that, as 
already mentioned, the Hamiltonian associated with the ground state solution (R, — i? 2 ,0) vanishes, one gets 

X 2 H = 2a 2 J (VR ■ Vo-i - R 3 <ri)£dZ 

+2b J (VR-Va 2 -R 3 a 2 )Zdt 
2V 2 f 

+ —J VR-Va 3 -R 3 a 3 )^ 



+ ^ J \AR\ 2 ^+^ J \V (R) 2 \ 2 + 

(50) 

Performing integration by parts in the integrals involving the oVs and using that R obeys AR — R + R 3, = 0, one gets 



\ 2 H = - J (2a 2 Ra x + 2bRa 2 + ^R<r 3 )i<lL 

" 2x2 V^ + y [\vim 



+1 J (|Ai?| 2 + i|V(i? 2 )| 2 )^ + --- (51) 

Using (|48|) , the first integral in the right-hand side of eq. (jSTj) is the difference Af = Af — J R 2 £,d£ between the actual 
plasmon number and its critical value for collapse when T = 0. Furthermore, a 2 = A 2 . One thus gets the following 
dynamical equation for the scaling factor A 



m 



_^^ A2+ A/--^), (52) 
i +m 2 A 2 V A 2 / K ' 



with mi = (1/2) / v 2 £d£ « 0.727, m 2 = (1/4) / £ 2 R 2 £d£ w 0.553 and m 3 = / (|Ai?| 2 + (l/2)|V(i? 2 )| 2 )^ w 10.785 
and Ji < 0. As in [3], the problem can be viewed as the motion of a particle in a potential that diverges positively 
as A — » and has a finite positive limit as A — > +oo. If T < A/" 2 / 4 1 | rn.3 , this potential has a negative minimum and 
A(i) oscillates between strictly positive minimal and maximal values. 

For an explicit numerical resolution, it is more convenient to consider y = A 2 that obeys 



m\ + m 2 y 



{ny 2 + Afy-Tm 3 ) (53) 



Under the condition T < A/" 2 /(4|7i|m3) , there exists positive y m and yu solutions of Tty 2 + Afy — r 777,3 such that y 
oscillates between y m = [Af - {Af 2 - A\H\Tm 3 ) 1/2 }/2\H\ and y M = [Af + {Af 2 - A\H\Ym z ) 1/2 ]/2\Hl consistent with 
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FIG. 2: Time variation of the scaling factor A (solid lines) for V = 5. 1CT 3 , H = -0.0430 and N = 0.240 (top) and V = 10" 3 , 
H — —0.0295 and J\f = 0.168 (bottom). For comparison, the corresponding evolutions in the case r = (dashed lines) are 
superimposed. 



the global existence demonstrated in Section III. Equation (|53[) is equivalently rewritten 

2 



Vtt = d y 



mi + TO22/ 



(Uy 2 +Ny~Tm 3 



For a given y(0), one should prescribe 



yt(o) 



7V(0)+A/V(0)-rm 3 
mi + m 2 y(0) 



1/2 



(54) 



(55) 



Figure [U shows the time evolution of the scaling factor A(i) for F = 5. 1(T 3 , H = -0.0430 and Af = 0.240 (top) 
and r = 10~ 3 , 7i = —0.0295 and Af — 0.168 (bottom), corresponding to the initial conditions Eq — 2.85 e~ r and 
Eq = 2.90 e~ r respectively, with uq — — \Eq\ 2 , Vq — 0. For comparison, we superimposed the corresponding evolution 
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FIG. 3: Profiles of the positive solution (ground state) of the equation A.^R^ — R^ + R^ 3 = and of the density 



when r = 0. In this case, A reaches zero in a finite time, with A 2 = N jm\ — Oq, consistent with the scaling A oc (t* —t) 
of the self-similar blowing-up solutions of the classical Zakharov equations. Note that, while near threshold the two- 
dimensional scalar model predicts the same leading- order profile for the pump wave as the nonlinear Schrodinger 
equation resulting from the subsonic approximation (slaved acoustic waves) , it does not lead to the same scaling law. 



V. THE TWO-DIMENSIONAL ELECTROSTATIC MODEL 

Let us now return to the electrostatic model, first in the case r = 0. From eqs. (Tr9"l) - (f2"U)) . we perform the 

rescaling = (1/A)S(£, r )e lMaA| £ |2/4 , n = (1/A 2 )A(£,t) and v = (1/A 2 )V(|, r). After neglecting as previously 
the r-derivatives of the rescaled functions, we get 

l£l 2 

AS-S + 6-^-S = 

4 

a 2 £(N) - AN = A|S| 2 . (57) 

The phase factors in eq. ([56")) introduce a serious difficulty in the sense that their expansions lead to an additional 
contribution in the perturbative calculation that is not necessarily associated with a well-posed problem, restricting 
de facto the present analysis to isotropic solutions for which the operator A _1 graddiv reduces to the identity and 
the phase factors cancel out. Denoting by S the radial component of S, one can then perform an analysis similar to 
that of the scalar model. 

The only difference between the isotropic electrostatic equations and the scalar model is the replacement of the 
scalar Laplacian in the equation for the electric field by the radial component of the vectorial Laplacian (which 
implies in particular that the electric field now vanishes at the center of symmetry). Denoting by the positive 
solution (fig. [3]) of 

A (i) i? (i)_ i? (i) +i? (i) 3 = (58) 

with ijW(O) = R^(oo) = (see [2(| for a review on this equation and its extensions that arise in the context of 
vortex solutions of the nonlinear Schrodinger equation), the rescaled density and velocity satisfy the same equation 
as in the scalar model, except that R is replaced by R^\ We thus recover eq. with coefficients now given by 

mi ~ 24.42, m2 w 8.14 and TO3 w 24.94. Furthermore J R^ £d£ w 7.69. As consequence, the scaling coefficient 
displays the same oscillatory behavior as in the scalar model. 
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VI. DIFFICULTIES IN THREE DIMENSIONS 



The three-dimensional problem displays specific difficulties even in the context of the scalar model. Indeed, although 
solutions of the classical scalar model (T = 0) with negative Hamiltonian appear to blow up in a self-similar way as 
in two dimensions, the perturbation analysis for T > is not straightforward in three dimensions. 

Proceeding as in two dimensions (but with different rescalings), we are looking for solutions of the form E(r, t) = 
A- 3 / 2 S(£,T)e^~ aAl/2 « 2 / 4 \ n(r,t) = \- 2 N(£,t), v(r,t) = \-^ 2 V^,t), n t (r,t) = X^^W^t), with £ = r/A with a 
profile depending only weakly on the rescaled time r. This leads to 

AS - S - NS + b^-S 
4 

_ r wAC 2 /4 e _ wA?2/4 = Q 
x z 

aC-W + ^W) - XAN = A\S\ 2 - J -A 2 N (60) 
Z A 

±d s (?d £ V) = -a(2N + ^N) = -W (61) 

+ (,d e V) - Xd^N = -d e \S\ 2 + —d^AN (62) 
Z A 

where a = -XtX 1 / 2 = -A r A~ 3/2 and b = (aX 1/2 ) T + (aX 1/2 ) 2 . From eqs. ((6TJ) and (J62J) , we have 

a 2 C{N) - XAN = A\S\ 2 - ^-A 2 N (63) 

A 

with now C = ^N^ + (13/2)£/V € + 7N. 

In the classical regime (r = 0), a has a finite (positive) limit ao as the collapse time is approached, corresponding 
to a scaling factor varying like (t+ — t) 2 ^ 3 . In this regime, the term XAN is negligible in eq. (|6"3"|) (supersonic regime) 
and S scales like a. However, as shown below, this ordering eventually breaks down in the presence of quantum 
effects. Indeed, as the collapse is arrested, At vanishes and so does a. This indicates that a systematic modulational 
theory analogous to what we developed in two dimensions, is not possible in three dimensions. In this context, we 
resort to limit our address of the problem to a semi-phenomenological approach based on the extension of asymptotic 
expansions outside their range of strict validity, with the hope to capture qualitative properties of the global dynamics. 

For this purpose, we assume that T is small enough for the quantum effects to start acting only after the system 
has reached the classical blowing up regime where a(t) is closed to its limit ao- In this asymptotic regime, one expects 
that there exists a period of time during which a remains sufficiently close to ao to allow a perturbative calculation. 
We are thus led to expand 

A r T a 2 — a 2 

S = a Q (Sa + bSi + -^S 2 + -tkS 3 + t-^t^ H n~ -S5) 

% X 2 Aoq ag 

N = N + bN 1 + \n 2 + ^N 3 + -^—Ni + °—^- N 5 

Qq A A(2q Qq 

V = a (V + bVi + \v 2 + ^-V 3 + -^V± + ?-^V 5 

Qq A AQq Qq 

+(a - a )V ). 
The lowest order terms (So, No, Vq) are solutions of 

AS - S a - N Q So = (64) 
C(N ) = A(S 2 ) (65) 
-C 2 ddt 2 Vo) = 2N Q + ^N . (66) 

^ + ^ = -^(l 5 o| 2 ) (67) 

Their profiles (not shown) are qualitatively very similar to those in two dimensions [l9l |. The corrections terms 
(Si,Ni,Vi) (i = 1,5) satisfy the systems 

ASi - Si - N Si - SoN, = F { (68) 

C(Ni) - 2A(S Si) = d (69) 

C%{?Vi) = Hi, (70) 
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with F 1 = -SfSo, Gi =0,F 2 = 0, G 2 = AN , F 3 = A 2 S , G 3 = 0, F 4 = 0, G 4 = -A 2 iV , F 5 = 0, G 5 = -C(N ), 
and Hi = —aa(2Ni + £d^Ni). We find that 5*5 = ^Sq and N5 = V5 = 0. The first and last term in the expansion of 

a 2 +a 2 

S thus combine giving 2a " Sq. Proceeding as in two dimensions, we substitute the expansions of S, N, V into the 
plasmon number and the ifamiltonian. We get 

a 2 r 

N = a a a 2 + 2 ai a 2 b + 2a 2 \ + 2a 3 T-^ + 2a 4 - + ■■■ . (71) 

where the coefficients cti = J SoSi^d^ are numerically estimated as ao « 1.99, ot\ ~ 1.94, a 2 ~ 0.35, a 3 « 51.36, 
a 4 w 34.44. On the other hand, 

X 2 H = I3 a 2 + (3 ia la 2 \ + f3 2 a 2 b + (3 3 X + p 4 T^ + [3 5 j + • • • (72) 



One first proves that 



i,-j{\VS \ 2 + n\So\ 2 + ^V 2 )e^ 0. (7:-!) 



For this purpose, we multiply eq. (|64[) by Sq and by £,{Sq)^ respectively, and integrate in space the resulting equations 
to get 

J (\VS \ 2 + N a S 2 + S 2 )ed£ = , (74) 
J (IVSol 2 + 3N S 2 + 3S 2 + ZN 0i S 2 )edt = . (75) 
Subtracting these two equalities and using (|66|) . we have 

[2S o 2 + y o d 5 (S o 2 )]£ 2 d? = 0. (76) 
Multiplying ([67)) by Vq and integrating in space, we get 

J v 2 e^ = - J (specie (77) 



Combining (JTBJ) and gives / S^ 2 d^ = \ fV£{, 2 d£, which after substitution in ([74]) leads to ([73]). 

Using eqs. (|64 ]) -([70 l) . the constants fa for i = 1, • • • 5 are given by /3\ = / ^S 2 1.17, /3 2 = / (x^o ~ 25 5i + 
Wi)£ 2 d£ « -1.17, /8s = JV2S S 2 + U V2 + ±N 2 )edt « 3.10, ft = J(-2S S 3 + TO^d? « -5.56, ft = 
/ (-25 5 4 + V V 4 + i|V7V | 2 K 2 dC « -2.81. 

We now eliminate the term proportional to a 2 b in H. using M and find 

a 2 

\ 2 TL = mia^a 2 \ + m 2 \ + m 3 T-^ + m 4 a 2 

A 

r 

+m 5 -+m 6 J\f. (78) 
A 

The new constants m, are defined as mi = ft w 1.17, m2 = ft — a^ft/ai ~ 3.31, 777,3 = At — d^ft/ai ~ 25.41, 
777 4 = — ao/?2/(2o;i) « 0.60, 777,5 = ft — ayft/ai « 17.95, mg = f3 2 /(2ai) ~ —0.30. 
Inserting that a — — AtA 1 / 2 , the effective ODE satisfied by A is 



, 2 _ 777,2 A 3 + m 3 Tal + m 5 TX + m 6 JVX 2 - HX 4 
X 2 (mia 2 l X 2 + m 4 A) 



K = X2/J -2 , I ~ 5 ( 79 ) 



In order to fix the parameter oq, we note that for r = 0, in the blowing up regime where A — > 0, one has 



a 2 = lim a 2 = lim A? A = -?**L ( 80 ) 

U A^0 A^O 777 4 
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0.25 




FIG. 4: Time variation of the scaling factor A for the three-dimensional phenomenological model, when Af — 5.64 and H = 
— 18.97, with F = 2. 10 -5 (top) and Af = 2.76 and Ti — —0.33, with the sale value of T as above. For comparison, the 
corresponding evolutions in the case r = (dashed lines) are superimposed. 



or using the definitions of m-4 and me, 



Af 



(81) 



Figure 4 (top) shows the evolution of the scaling factor X(t) for Af = 5.64 and Tt = —18.97 (corresponding to 

an initial conditions Eq = ce~ r , uq = — |i?o| 2 , and vq = with c = 6, for the original Zakharov equations) and 
r = 2. 10~ 5 . For comparison, the evolution in the absence of quantum effects is also displayed. As in two dimensions, 
we observe that quantum effects induce a periodic behavior. In simulations with smaller T, the minimum of A is, as 
expected, getting smaller. When keeping the same value of T, one uses H = —0.32 with Af — 2.76 by taking c = 4.2, 
the maximum of A is slightly increased, while the period of the oscillation gets significantly longer (Fig. 4 bottom), 
but the global behavior remains very similar. Note however that the dynamics is much faster than in two dimensions. 
This effect is already visible on the singularity time at T = 0, a regime for which eq. (f79")) is asymptotically exact. 
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This suggests that the parameters we used in three dimensions (even for the bottom panel of Fig. 4) correspond 
to a regime significantly distant from the threshold conditions that in this case are not as easily characterized as 
in two dimensions. As stressed at the beginning of this section, the present description of the three-dimensional 
problem is however to be viewed as heuristic, as the asymptotics clearly breaks down before the scaling factor X(t) 
reaches its minimum. It nevertheless predicts a behavior consistent with the arrest of collapse (the proof given in 
Section I is easily transposed to the scalar model). It also shows a periodic dynamics whose origin is expected to be 
generic. It indeed results from a competition between wave focusing that occurs when A is not yet small enough for 
the quantum effects to act efficiently, and the subsequent evolution that takes place when the influence of the latter 
perturbations dominates the dynamics and leads to defocusing until the moment where, A becoming large enough, 
their influence becomes again subdominant, thus permitting an efficient self- focusing. Validation of the model would 
require comparisons with direct simulations of the scalar model, an issue that is outside the scope of the present 
paper. Compared with the Rayleigh-Ritz method, the present approach should provide a better description of the 
solution profile. It also more clearly points out the conditions of applicability of modulation methods. 



VII. CONCLUSION 



The influence of quantum effects on the Langmuir wave dynamics provides an interesting example of the action 
of an additional dispersive effect on the phenomenon of wave collapse, although in realistic situations the Zakharov 
description is supposed to break down before quantum effects become relevant. Arrest of collapse was predicted Q 
in the adiabatic regime where the density is slaved to the wave amplitude, using a Rayleigh Ritz method. Here, the 
result is rigorously established for the full quantum Zakharov equations, by combining the conservation of the plasmon 
number and of the Hamiltonian with estimates based on a Gagliardo-Nirenberg inequality. These invariances are also 
used to develop a systematic perturbative expansion in order to capture the influence of weak quantum effects for 
initial conditions slightly above the singularity threshold for the classical problem. Restricted for technical reasons to 
isotropic solutions (at least in the focusing region) , the analysis is carried out in two space dimensions corresponding to 
the critical dimension for collapse of the cubic nonlinear Schrodinger equation. The difficulty of extending the analysis 
to three dimensions points out the importance of the proximity of criticality to satisfy the delicate balances involved 
in a systematic asymptotic theory. We thus resorted in this case to develop a semi-phenomenological approach. 
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